A Model for Quantum Jumps in Magnetic Resonance Force Microscopy 
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We propose a simple model which describes the statistical properties of quantum jumps in a single- 
spin measurement using the oscillating cantilever-driven adiabatic reversals technique in magnetic 
resonance force microscopy. Our computer simulations based on this model predict the average time 
interval between two consecutive quantum jumps and the correlation time to be proportional to the 
characteristic time of the magnetic noise and inversely proportional to the square of the magnetic 
noise amplitude. 
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I. INTRODUCTION 

Recent achievements in magnetic resonance force mi- 
croscopy (MRFM) promise a single spin detection in the 
near future 0. One may expect that single spin signal 
will represent a random sequence of quantum jumps. The 
important problem for the theory is modeling of quantum 
jumps in MRFM and the computation of their statistical 
characteristics. In this paper we propose a simple model 
which describes quantum jumps in MRFM single spin 
detection. We consider the oscillating cantilever-driven 
adiabatic reversals technique (OSCAR) which currently 
is the most promising approach for single spin detection 
[3. In OSCAR the cantilever tip (CT) with a ferromag- 
netic particle oscillates, causing the periodic adiabatic 
reversals of the effective magnetic field on spin. The 
spin follows the effective magnetic field causing a tiny 
frequency shift of the CT vibrations which is measured. 

In the second section we consider the Schrodinger dy- 
namics of the CT-spin system which underlines our model 
of quantum jumps. In the third section we describe our 
model and present the results of the computer simula- 
tions. In Conclusion we discuss our results. 



II. SCHRODINGER DYNAMICS OF THE 
CT-SPIN SYSTEM 

We consider a vertical cantilever with a ferromagnetic 
particle attached to the CT and oscillating along the x- 
axis which is parallel to the surface of the sample. (See 
Fig. The Hamiltonian of the CT-spin system in the 
system of coordinates which rotates with a circularly po- 
larized r/ field can be written as 



+ xl) + eS,.+2r^XcS, + A{t)S,. (1) 



Here pc is the effective momentum of the CT, Xc is its 
coordinate, the first term describes the effective energy 
of the CT; the second term describes the interaction be- 
tween the spin and the r/ field; the third term describes 
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FIG. 1: Single spin OSCAR MRFM setup. Bext is the ex- 
ternal permanent magnetic field; Bi is the rotating rf field; 
rh is the magnetic moment of the ferromagnetic particle; s is 
a spin near the surface of the sample. 



the interaction between the CT and the spin; and the last 
term describes the interaction between the spin and the 
random magnetic field which causes a deviation of spin 
from the effective magnetic field when the latter passes 
through the transverse x — y-plane 0, Q ■ AH quantities 
in (1) are dimensionless: Pc and expressed in terms 

of the "quantum units" Aq and Pq 



Ao - {huJk.Y'^, Po - h/Xo 



(2) 



The dimensionless parameters e, 77, and A(t) are defined 
as 

e = 7EiK, rj = -f{h/k,LU,y/^\dB,{x)/dx\/2, (3) 



A(r) ^jAB,{t)/uJc, T^u^t. 

Here is the z-component of the "regular mag- 
netic field" which includes Bf,xt and the regular dipole- 
magnetic field produced by the ferromagnetic particle at 



the location of the spin; /^B^^ir) is the z-component of 
the random field with zero average value; 7 is the mag- 
nitude of the spin gyromagnetic ratio: kc and ojc are the 
effective spring constant and fundamental frequency of 
the CT; and r is the dimensionless time. 
Using the parameters presented in Q 

uJc/2tt = 6.6kHz, fc^ = 6 x IQ-^N/m, (4) 

Bi = O.SmT, \dB,Xx)/dx\ « 4.3 x lO^T/m, 

we can estimate all parameters in (1) except of A(r): 

Xa = 85/to, Pa = 1.2 x IQ-^^Ns, 77 = 0.078, s = 1270. 

(5) 

To simplify computer simulations wc considered the func- 
tion A(t) to be a random telegraph signal with two 
values ±A. The time interval between two consecutive 
"kicks" of A(r) was taken randomly from the interval 
(tq — St, To -f- St) , with the average time interval, tq , close 
to the Rabi period t/j 

TR = ujc{2n/jBi) = 27r/£ = 4.95 x 10"^ (6) 

We choose the initial state of the CT to be a coherent 
quasiclassical state, and the average spin to be pointed 
along the "regular" effective magnetic field with the com- 
ponents {e,0,2r]Xc). 

Below we describe briefly the results of our computer 
simulations. Our simulations reveal the formation of a 
Schrodingcr cat state for the CT: the probability distri- 
bution function P{x,t) = '^'^ [x , t)"^ {x , t) splits into two 
peaks. Similar to our previous computations ^ IS 
one peak corresponds to the average spin pointing in the 
direction of the effective magnetic field, while the other 
one corresponds to the opposite direction of the average 
spin. The two peaks oscillate with slightly different pe- 
riods due to the back action of the spin as expected in 
the OSCAR technique [tI (sj . Unlike our previous papers 
^ in IS I the appearance of two peaks is not connected 
with the initial deviation of the spin from the direction of 
the effective field. It is induced by the action of the ran- 
dom field A(r) which causes this deviation in the process 
of the spin-CT dynamics. It was shown in 0,1^1 that the 
interaction with the environment quickly destroys quan- 
tum coherence between the two peaks: the Schrodinger 
cat state of the CT quickly transforms into a statistical 
mixture of the two possible CT trajectories. Thus, the 
random magnetic field A(r) causes the quantum jumps: 
spin may "jump" to the opposite direction relative to 
the effective magnetic field and, correspondingly, the CT 
may slightly change the period of its oscillations. 

Unfortunately, our present simulations consume too 
much computer time to be able to reveal more than 
one spin jump. This approach, clearly, does allow one 
not to study statistical characteristics of quantum jumps. 
Consequently we developed a simplified model which de- 
scribes statistical properties of the spin jumps. 




jump jump 

FIG. 2: (a) Distribution function of time intervals between 
two consecutive quantum jumps for A = 100 tq = 0.01 and 
10^ kicks; the solid line is a fit with Td = 32; (b) Enlargement 
of (a) 



III. SIMPLE MODEL OF QUANTUM JUMPS 

We assume that every "kick" provided by the function 
A(r) is followed by the collapse of the wave function. 
Before the kick, the spin points in (or opposite to) the 
direction of the effective magnetic field. After the kick 
there appears the finite angle (AO) between the new di- 
rection of the effective field and the average spin. The 
probability for the spin to "accept" the "before-kick" di- 
rection relative to the new effective field is cos^(A0/2). 
The probability to "accept" the opposite direction, i.e. 
the probability of a quantum jump is sin^(A0/2). (It 
is clear that a significant probability of a quantum jump 
appears only when the effective field passes the transver- 
sal X — y-plane. In the transversal plane, the effective 
field has the minimum value. Thus, after every kick of 
the random field our computer code decides the "fate" 
of the spin in accordance with the probabilities of two 
events: to restore the previous direction relatively to the 
effective field, or to experience a quantum jump. In our 
model the CT experiences harmonic oscillations 

2^1^' (''') = cos(l ± 5uj)t, (7) 

where (±) correspond to two CT trajectories with the 
spin pointing in the direction of (or opposite to) the cor- 
responding effective field, and Slo is the CT frequency 
shift. The components of the effective field are given by 

(£,0,277Xe(T) + A(r)). (8) 

From the experimental data 0] for the CT amplitude 
X„i = 10 nm we obtain Xm = 1-2 x 10^. The frequency 
shift, Suj, can be estimated as Q 

Slo^Alu/uJc^ "^^^l =4.2x10-^ (9) 
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FIG. 3: Dependence of the average time intervals between 
two consecutive quantum jumps on ro/A'^. The fuU line cor- 

, the squares repre- 
4.2 X 10"'', the cir- 



responds to the value Xm ~ 1.2 x 10 
sent Slo — 4.2 x 10~*, the crosses Sui 
cles Suj = 4.2 X 10" 
value Xm = 7.2 x 10 
ing parameters A and to in the ranges 
0.001 < TO < 1. 



The dashed line corresponds to the 
^. Data have been obtained by vary- 
10 < A < 300 and 



Note that our model contains two important simplifi- 
cations: first, we assume that the wave function collapse 
occurs immediately after the "kick" of the random field. 
Thus, we ignore the finite time when the spin-CT system 
is in an entangled state. Second, in a real situation the 
deviation of the spin from the effective field is a "quasi- 
resonance" process caused by the cantilever modes whose 
frequencies are close to the Rabi frequency. In our model 
this deviation appears as a result of the "kick" of the 
random field. 

Below we describe the results of our simulations. 
Fig. 121 demonstrates a typical probability distribution 
of time intervals, Tjump, between two consecutive quan- 
tum jumps. The probability distribution is a sequence of 
sharp peaks at Tjump = Tn = mr with the Poisson-like 
amplitude 



-P(t„) ^ exp(-r„/r<j). 



(10) 



(Certainly, P(tj, 



at r < To — St.) The sharp 



peak appears as the probability of the quantum jump is 
significant when the spin passes through the transversal 
plane, i.e. every half-period of the CT oscillation which is 
equal to tt. The average value of the time interval (r- 
was found to be 



jump / 



i^jump) 



(11) 



with an error less than 3%. The standard deviation is 
equal to Td with the same accuracy 



{{'''jump) {Tjump) ) ^ ~ Td- 



(12) 



We studied the dependence of the average value 
Tjump) on the parameters of our model. We have found 
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FIG. 4: Correlation function of the CT frequency shift C(r) 



A = 50, TO = 10~^ triangles A = 50, ro = lO'^ 



In all 



cases 5r — ro/4. The dashed curves show the exponential 
approximation of the correlation function exp(— t/tc) with 
Tc = 23.91, 95.81, 1179.15, respectively. 



that {Tjump) does not depend on St or Slu. (We varied 
St from to tq and changed Slo up to one order of mag- 
nitude.) At a fixed value of the amplitude x„i the value 
of (Tjump) is approximately proportional to tq/A^. Fig. 
13 demonstrates this dependence. 

The best fit for the numerical points in Fig. |3|is given 

by 



ln(r. 



jump / 



:p + gln(To/A2). 



(13) 



For Xm = 1-2 x 10^ we have p = 17.9, q = 0.993. For the 
six fold value Xm = 7.2 x 10^ we obtained the same value 
of q and p — 19.743. If we estimate the amplitude of 
the random CT vibrations near the Rabi frequency as 1 
pm, then A = 1.8. Putting tq = tr and the experimental 
value for tr (6) we obtain LUc{Tjump) — 2.3 s for Xm — 
lOnm and uJc{Tjump) = 14.5 s for x,„ = 60nm. These 
values are close to the experimental characteristic times 
of 3 s and 20 s reported in . 

Next we computed the correlation function for the CT 
frequency shift 



C{t) 



{{Slo (t) - Suj) (Suj (t - 



Slo)) 



({Slo (t))^) - Slo 



(14) 



where Slo — (SLo{t)) — 0, and (...) indicates an average 
over time. 

In Fig. 4 we show the correlation function C(r) for three 
different values of parameters as indicated in the legend. 
As one can see, the general behavior is well described 
by the exponential function (indicated by dashed lines in 
Fig. 4) exp(— t/tc). The relation between the correlation 
time Tc and (Tjump) was found to be {Tjump) — 2.5rc. 
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IV. CONCLUSION 

We developed a simple model which describes quan- 
tum jumps of spins in the OSCAR MRFM technique and 
which allowed us to compute the statistical characteris- 
tics of quantum jumps. In our model the average time 
interval {rjump) between two consecutive jumps and the 
correlation time Tc are proportional to the characteristic 
time of the magnetic noise, and inversely proportional to 
the square of the magnetic noise amplitude. Using exper- 
imental parameters [l| and a reasonable value for the CT 
noisy amplitude we obtained the value of {rjump) which 



is close to the experimental value of the characteristic 
time of OSCAR MRFM signal. 
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